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Abstract 

We study the linear stability with respect to lateral perturbations of free surface films of polymer 
mixtures on solid substrates. The study focuses on the stability properties of the stratified and 
homogeneous steady film states studied in Part I [U. Thiele, S. Madruga and L. Frastia, Phys. 
Fluids 19, 122106 (2007)]. To this aim, the linearized bulk equations and boundary equations are 
solved using continuation techniques for several different cases of energetic bias at the surfaces, 
corresponding to linear and quadratic solutal Marangoni effects. For purely diffusive transport, an 
increase in film thickness either exponentially decreases the lateral instability or entirely stabilizes 
the film. Including convective transport leads to a further destabilization as compared to the 
purely diffusive case. In some cases the inclusion of convective transport and the related widening 
of the range of available film configurations (it is then able to change its surface profile) change the 
stability behavior qualitatively. We furthermore present results regarding the dependence of the 
instability on several other parameters, namely, the Reynolds number, the Surface tension number 
and the ratio of the typical velocities of convective and diffusive transport. 
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I. INTRODUCTION 



The main stages of structure formation and dewetting in thin films of simple liquids are 
generally well understood^-. Films of polymers rupture forming holes, and evolve towards 
a network of liquid rims that may decay subsequently into small droplets^. Moreover, 
using physically /chemically heterogeneous substrates it is possible to create ordered soft 
matter films^— . For single component thin films there are two main pathways of dewetting: 
spinodal surface instability and heterogeneous nucleation, both driven by effective molecular 
interactions between the substrate and the film surface like, e.g., van der Waals forces. Both 
mechanisms may result in the rupture of the film leading to various patterns and have been 
studied extensively in numerical and experimental works.— _ — 

The dewetting of thin films composed of binary mixtures is more involved. Experiments 
on mixtures have reported dewetting mechanisms that differ from the ones for simple liquids. 
In particular, Yerushalmi-Rozen et a/.— discuss the phase-separation induced dewetting of a 
polymer blend consisting of deuterated oligomeric styrene and oligomeric ethylene-propylene. 
They assume that first there occurs a separation of the two components of the blend. The 
phase separation is followed by a dewetting process that is characterized by the formation 
of holes at a dewetting front moving inwards from the sample edge. This route of dewetting 
differs from the classical ones for single component films by its short induction time (at 
comparable film thickness), the properties of the front, and its morphological characteristics. 

The evolving gradients of concentration of the two components of the binary mixture along 
the film surface originate surface tension gradients that are responsible for an additional 
flow in the vicinity of the decomposition/dewetting front. It is thought that this promotes 
the acceleration of the formation of holes. However, no detailed theoretical description 
is available at present. The relevance of interfacial phenomena in binary mixtures has as 
well been reported in spin-cast polymer blend films, where interfacial instabilities lead to a 
horizontal phase separation^. 

Another element relevant to the dynamics of films of binary mixtures is the energetical 
influence of the surfaces. In particular, an energetically biased surface is rapidly enriched in 
the preferred component and may become the 'source' of so-called spinodal decomposition 
waves^ - — . Jones et alM, for instance, study the spinodal decomposition of critical mixtures 
of poly(ethylene-propylene) and its per-deuterated variant in presence of a substrate with 
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a preferential attraction for the poly(ethylene-propylene). They find that that composition 
waves origin at the surfaces and propagate into the bulk. Their wave-vectors are oriented 
normal to the surface. 

The interplay between phase separation, surface tension gradients, and surface properties 
make the dynamics of films of mixtures extremely rich and its theoretical description chal- 
lenging. The dynamics of a fluid binary mixture is often described by the so called model-H, 
which couples transport of the mass of one component (convective Cahn-Hilliard equation) 
and momentum (Navier-Stokes-Korteweg equations)2§*2Z. It has been extensively studied 
for various isothermal and non-isothermal settings^— . None, however, involves beside the 
diffuse internal interface a free liquid-gas interface, i.e., a free surface. This implies that 
existing theories can not be readily applied to thin film experiments with polymer blends 
whose evolution is governed by the driving influence of the dynamics of and at the free 
surface. 

In an alternative approach Clarke^ 1 ^ constructs and investigates a simple thin film model 
for the coupled evolution of film thickness profile and mean concentration. It is, however, 
based on the assumption that the films show no vertical concentration profile. This implies 
that it can not be applied to the case of the lateral stability of layered films we are interested 
in here. 

The present series of works aims at the development and application of a model-H for 
systems involving free surfaces. Note that the incorporation of hydrodynamic flow is nec- 
essary even in the case of extremely slow cree ping flow. Otherwise no evolving surface 
deflections can be described. The first part (Ref. 1371 ) discusses the inclusion of free evolving 
surfaces into model-H and determines basic stratified film states for various types of ener- 
getic bias at the surfaces. In particular, the framework of phenomenological non-equilibrium 
thermodynamics is used to derive a generalized model-H coupling transport equations for 
momentum, density and entropy that is then simplified for an isothermal setting, vanishing 
interface viscosity and simplified internal energies. A comparison with literature results and 
a variational derivation of the static limiting case clarify the issue of defining pressure and 
chemical potential. Furthermore, boundary conditions at the solid substrate and the free 
interface are introduced. After non-dimensionalisation it is shown that the dimensionless 
numbers entering the boundary conditions for the Cahn-Hilliard and the Korteweg-Navier- 
Stokes equations are closely related and can not by any means be chosen independently of 
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each other. Physically this means that the energetic bias with respect to decomposition 
at a free surface is intrinsically coupled to a solutal Marangoni effect. Neither of the two 



effects can be considered independently of the other one. Ref. l37| continues with an analysis 
of steady base state solutions for laterally homogeneous films of decomposing mixtures. A 
plethora of stratified solutions is determined and ordered for various types of energetic bias 
at the surfaces using continuation techniques and symmetry arguments. 

This present second paper presents the analysis of the lateral stability of the laterally 



homogeneous steady films states obtained in Ref. l37|, i.e., it determines the stability of 
homogeneous and layered films with respect to modulations in thickness and/or composition 
along the substrate. A similar approach (however, with other boundary conditions than used 
here) was recently employed to obtain instability thresholds for long-wave instabilities^. 
The paper is organized as follows. In Section [Til we review the non-dimensional model-H 
including appropriate boundary conditions at the solid substrate and the free film surface. 
In Section II III equations are linearized for base states corresponding to homogeneous and 
stratified films. We as well discuss a realistic range of parameters for polymer blends to be 
used in the calculations. The stability results for homogeneous films are presented in Section 
IIVI for neutral and energetically biased surfaces. The stability of stratified films for neutral 
and biased surfaces is presented in Section |VJ Finally, conclusions are given in Section IVIl 

II. DIMENSIONLESS MODEL 

We consider a film of a binary mixture on a horizontal homogeneous solid substrate. 
The film has a free surface that may evolve in time. The system is two-dimensional and 
infinitely extended in the horizontal direction. The origin of the Cartesian frame is fixed to 
the substrate. The dimensionless governing equations expressing the balance of mass and 



momentum, and the continuity equations are re-derived and discussed in Ref. |37|. They read 
Ps 



<9v „ 
-+V-W 



-V.{(Vc)(Vc)+p c ,I} + ^Av, (1) 

Re 



d tC + v • Vc = -V • {V[Ac - dj(c)}} , (2) 

and 

V ■ v = 0, (3) 
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respectively, with the operators V = (d x ,d z ), A = (d 2 + d 2 ), and the velocity field v = 
(u,w). The composition field c = c\ — c<i represents the difference of concentration of the 
two components of the mixture, p c g = p — (c + l)Ac — (Ac) 2 /2 is an effective pressure 
comprising all diagonal terms of the stress tensor, and p is the 'normal' pressure. The local 
bulk free energy is assumed to correspond to the simple quartic potential f(c) = (c 2 — l) 2 /4. 



To make the equations ([T]) and ([2]) dimensionless, in Ref. |37| the scales I = C y o c j E, U = 
M E/IC 2 , t = l/U = I 2 C 2 / (ME) and P = E are used for length, velocity, time and pressure 
respectively. The length / corresponds to the thickness of the diffuse interface between the 
two phases of the mixture and is determined from the coefficient a c of the gradient term of 



the convective Cahn-Hilliard equation (cf. Ref.l37l). the energy scale E and the concentration 



C at the binodal. The parameter M is the diffusion coefficient in the Cahn-Hilliard equation. 
With these scales, two dimensionless numbers appear in the bulk equations (JT]) and ©, the 
pressure number Ps = p M 2 E 2 / C 6 cr c and the Reynolds number Re = M Ep/rjC 2 . Note that, 
as the diffusion constant is D = ME and U ~ D/l, the used Reynolds number can be seen 
as an inverse Schmidt number Sc= r]/pD. The ratio Re/Ps=cr c C 4 /rjME will turn out to be 
the most important bulk parameter for the extremely slow creeping flow we are interested 
in (see below Section [HI Bj) . It corresponds to the ratio of the typical velocity U — cr/rj of 
the viscose flow driven by the internal 'diffuse interface tension' o = a c /l and the typical 
velocity of diffusive processes U ~ M E/l = D/l. This implies that Re/Ps can as well be 



seen as a Peclet number Pe=U'l/D or as an inverse Capil 
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and 



ary number Ca 1 = a/Urj. Such 



31 



alternative scalings are used, for instance, in Refs. 

Equations (pD) and (j2j) are supplemented with boundary conditions for the concentration 
and velocity fields.— We first present the conditions for the concentration field. At the solid 
substrate (z — 0) one prevents a diffusive mass flux through the substrate and allows for an 
energetic bias, i.e., 

d z [(d xx + d zz )c-d c f(c)] = (4) 
-d z c + Sd c f-(c) = (5) 

respectively. The energy bias is f~(c) = 7 S + a~ c + b~c 2 /2, where S7 S is the dimensionless 
solid-liquid interface tension at c = 0. Note that it does not influence later calculations 
as only d c f~(c) enters. Parameters a~ and b~ model preferential adsorption of one of the 
species at the substrate and changes in the mixing/ demixing behaviour of the species at the 



substrate, respectively. The dimensionless parameter S= Jq/IE is the dimensionless surface 
tension of the liquid-gas interface, and 70 is the reference surface tension at c = 0. Similar 
conditions are applied at the (curved) free surface (z = h(x,y,t)), i.e. 

[d z - (d x h)d x ] [(d xx + d zz )c - d c f(c)} = (6) 
[d z - (d x h)d x ] c + Sd c f + (c) yj 1 + (d x h) 2 = (7) 

with f + (c) = 1 + a + c + b + c 2 /2. Parameters a + and b + quantify preferential adsorption of 
one of the species and changes in the the mixing/demixing behaviour at the free interface, 
respectively. 

The boundary conditions for the velocity field are the no-slip and no-penetration condition 
at the substrate, i.e., v = w = at z = 0. At the free surface the conditions result from the 
balance of tangential forces 



- [d x c + (d x h)d z c\ [d z c - (d x h)d x c] (8) 
Ps 

+ Re ^ Uz + Wx ^ 1 ~ h ^ + 2 ( Wz ~ Ux ) hx \ 
= $y/l + K[d x + {d x h)d z ] f + (c) 



and normal forces 



[d z c - (d x h)d x cf - p cS (9) 



1 + hl 
Ps 2 

+ \ u x h l + w z~ h x (u z + w x j] 

Re 1 + hi L 



Sf + (c)d x 



d T h 



L(i + ^) 1/2 . 

Notice that the tangential gradient of f + (c) enters the tangential force balance corresponding 
to a solutal Marangoni force. Furthermore, at the free surface one has a kinematic condition 
assuring that the free surface follows the velocity field 

d t h = w — ud x h. (10) 



37 embarks on an 



After deriving the system of equations and boundary conditions, Ref. 
extensive study of quiescent, vertically homogeneous and vertically stratified base state 
solutions that are all homogeneous (or translationally invariant) with respect to the direction 
parallel to the substrate. The corresponding concentration profiles cq(z) are obtained by 
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solving the steady bulk equation (steady version of Eq. (j2J)) 



= d zz [d zz c - d c f(c)\ C0 ] 



(11) 



with the boundary conditions 



= d z [d zz c - d c f(c)\ CQ ] 



at 



z = 0, ho 



(12) 



and 



= ±d z c + Sd c f ± (c)\ Co at z = 0, h . 



(13) 



Here we are concerned with the stability of the base state concentration profiles cq(z) with 
respect to perturbations along the direction parallel to the substrate. The rational behind 
this approach is the experimental observation that thin films of a decomposing mixture first 
stratify vertically on a relatively short time scale and then develop on a slower time scale a 
horizontal structure. Here, we want to capture the characteristics of the latter process by 
the linear stability analysis. 

III. LINEARIZED EQUATIONS FOR VERTICALLY STRATIFIED AND HOMO- 
GENEOUS FILMS 

1. General ansatz 

To analyze the stability of the quiescent base states with respect to infinitesimally small 
perturbations, we write the general solution of the problem in the form v = vo + £Vi, 
Pes = Po + £pi, c = c + ec 1 , and h = h + shi, with v = and p = — (<9 z c ) 2 . The 
fields evi, epi, ed\, and eh\ denote the infinitesimal perturbations of velocity, pressure, 
concentration, and thickness fields, respectively. The small parameter e will be used to 
order terms in the series expansion. The perturbations are decomposed into a sum of normal 
modes (vx,px, ci, hi) = (vi(z),pi(z), c±(z), hi) exp((3t + ikx), where (3 is the growth rate and 
k the lateral wavenumber. Using this ansatz in Eqs. (CQ) to (TlPl) we obtain the linearized 
convective Cahn-Hilliard equation 



(3ci + wid z c 



{dzz ~ k 2 ) [(d zz - k 2 )d - cid cc f\ co ] . 



(14) 
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and the linearized momentum equation 

/3Ps Vl = -V-{[(Vc )(Vc 1 ) + (Vc 1 )(Vc )]+p 1 l} 

+^(-k 2 + d zz )^ (15) 
Re 

where the tensor [(Vco)(Vci) + (Vci)(Vc )] = [(0, ikcid z co), (ikcid z co, 2(d z ci)(d z co))], i.e. 
splitting the velocity field in its components we obtain 

/3Psui = —ikd z (cid z co) — ikpi 

+^ e (-k 2 + d zz ) Ul (16) 

f3P SWl = -[-k 2 c 1 d z c + 2d z ((d z c )(d z c 1 ))}-d zPl 

+^(-k 2 + d zz ) Wl . (17) 
Re 

The incompressibility condition leads to 

%ku\ + d z w\ = 0. (18) 

The linearized boundary conditions assuring zero mass flux through the substrate at z = 
and through the free surface at z = h are both of the form 

= d z [(-k 2 + d zz ) Cl - Cl d cc f\ CQ ] . (19) 

The linearized energy bias conditions for the concentration are 

-d zCl + S Cl d cc f-\ C0 = at z = (20) 
and d z Ci + Sci<9 cc / + | co = at z = h . (21) 

The non-slip condition for the velocity at the substrate becomes 

Ul = w 1 = 0. (22) 

and the linearized tangential and normal force balance at the free interface are 

{d z u x + ik Wl ) = 0, (23) 

Ke 

2(d z c )(d z c 1 )+ Pl -2^d z w 1 = Sf + \ Co k 2 h, (24) 
respectively. Finally, the kinematic condition reads 

w x = (3hx. (25) 
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Note that the r.h.s. of Eq. (|23|) is exactly zero following from the boundary condition Eq. ({TBI) 
of the base state. Physically this means that the Korteweg and the Marangoni stress at the 
free surface exactly compensate, i.e. all the tangential driving force is already contained in 
the linearized bulk equation. The situation would be different if we were to allow for a 
dynamic surface tension different from the used static one. 



2. Eigenvalue problem for stratified film 

To carry out the stability analysis of the linearized model-H we eliminate the pressure in 
the momentum equation and boundary conditions, and write the linearized model-H as an 
eigenvalue problem of the form 

- [(k 4 - 2k 2 d zz )c l - (d zz - k 2 )( Cl d cc f\J] , (26) 



and 



Re 

dzzzzWi = (3 Re {d zz - k 2 )w! + — k 4 c x d z CQ 

Jr S 

Re 

- — 2k 2 d z [(d z c )(d zCl )} 
Re 

+ k'dzzicAco) + 2k 2 d zz w 1 - k 4 w u (27) 

Jr S 



with the boundary conditions 



Wl = d z wi = at z = 0; (28) 

Ps 

— (d zz Wi + fc 2 u>i) = at z = h ; (29) 
Re 



and 



k 2 (dzC )(dz Cl ) - k 2 Cl d 2 z c + ^(d zz -3k 2 )d zWl 

Re 

-pPsd z W! = S/ + | C0 £; 4 /ii at z = h ; (30) 
where h\ = Wi//3. Furthermore, 

d z [(-k 2 + d zz )ci - cid cc f\ CQ ] =0 at z = 0, h ] (31) 

- d z cx + Sci9 cc /"| C0 = at z = 0; (32) 
<9 2 ci + Sci9 cc / + | Co = at z = h . (33) 
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3. Eigenvalue problem for homogeneous films 



For homogeneous films the concentration does not depend on z, i.e., the base state is 
uniform along both, the x— and z-direction. In consequence, the eigenvalue problem reduces 
to the bulk equations 

d zzzzCl = -/3 Cl - (/c 4 - 2k 2 d zz ) Cl - (d zz - k 2 )( Cl d c J\J, (34) 

and 

d zzzz w 1 = (3Re (d zz - k 2 )w 1 + 2k 2 d zz w 1 - fc 4 Wi, (35) 
with the boundary conditions 

w\ = d z wi = at z = 0; (36) 

Ps 

^(d zz Wi + k 2 wi) = at z = h ] (37) 

Ps Ps 

— (d zz - k 2 )d z wi - (3Psd z W! - — 2k 2 d z w x 
Ke tie 

— Sf + \ Co kwx/f3 at z = ho] (38) 
d z [{-k 2 + d zz ) Cl - Cl d cc f\ C0 ] = at z = 0,h ] (39) 

- d z c x + Sci9 cc /"| C0 = at z = 0] (40) 
and d z c\ + Sci<9 cc / + | Co = at z = h . (41) 

Inspection of the equations shows that the perturbations in the concentration and in the 
velocities do decouple entirely, i.e., the linear stability problem for the concentration field in 
a homogeneous film reduces to the one resulting from the Cahn-Hilliard equation in a slab^. 
The decoupled velocity perturbations are all damped out. 

A. Numerical technique 

In the present work the calculations of the base states and their linear stability are carried 
out using numerical continuation techniques bundled in the package AUTOi^r— . 

Continuation techniques allow to obtain solutions of a problem for a given set of control 
parameters by 'extrapolation' from known solutions that are nearby in the parameter space. 
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In particular, the set of equations for the steady state (fTTl) and the for real linear pertur- 
bations (]26p -(l27 p can be written as a 10- dimensional dynamical system y'(z) = f(y(z),A), 
where y = (co, c 0z , c u c lz , c lzz , c lzzz , w u w lz , w lzz , w lzzz ) and A denotes a set of control pa- 
rameters (in our case Re, Ps/Re, S, /3, k, c). This system of ordinary differential equations 
together with the boundary conditions at z = and z = h ( I28|) - fl33|) and one integral condi- 
tion (mass conservation) is discretized in space and the resulting algebraic system is solved 
iteratively. The package AUTO^ uses the method of orthogonal collocation to discretize 
solutions, where the solution is approximated by piecewise polynomials with two collocation 
points per mesh interval. The mesh is adaptive to equidistribute the discretization error. 
Starting from known solutions, AUTO tries to find nearby solutions to the discretized sys- 
tem by using a combination of Newton and Chord iterative methods. Once the solution has 
converged, AUTO proceeds along the solution branch by a small step in the parameter space 
defined by the free continuation parameters and restarts the iteration. Boundary conditions 
and integral conditions require additional free parameters which are determined simultane- 
ously and are part of the solution of the differential equation. The package AUTO is limited 
to the continuation of ordinary differential equations, and has been successfully applied to 
other thin film problems like, e.g., the determination of dispersion relations for transverse in- 
stabilities of advancing liquid fronts and ridges^; the determination of steady and stationary 
thickness profiles and their linear stability for sliding^— , running^ or depinning^ droplets; 
the analysis of steady film profiles for epitaxial growth^; and the analysis of surface waves 
in falling film problems^— . 

The main difficulty is usually to provide a starting solution for the continuations. Here, 
we follow different sequences of steps starting from trivial analytically known solutions. 
For n = branches (for definition see below) we (i) increase the thickness starting from a 
homogeneous solution of small thickness, (ii) increase the energetic bias at the surfaces if 
needed, (iii) look for branching points when increasing the linear growth rate for some fixed 
wave number, (iv) calculate the dispersion relation and its maximum, (v) and follow the 
maximum in a multi-parameter continuation. For b ^ branches we (i) start with some 
fixed thickness, increase the energetic bias at the surface if needed and identify branching 
points where non-trivial solutions emerge (ii) increase the thickness following the various 
emerging branches, (iii)-(v) as above. In this way the linear stability of the various stratified 
films can be determined in a rather effective manner. 
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B. Parameter range 



Due to the large number of dimensionless parameters needed to describe a film of a bi- 
nary mixture with a free surface, we proceed to discuss the stability results in a restricted 
range of parameters that is of particular experimental interest. For instance, the polymer 
blend Polystyrene/Polyvinylmethylether (PS/PVME) has been widely used in thin film ex- 
periments, without^ 1 ^ and with^ 3 - the presence of an external electric field. We will use this 
mixture as a reference for our calculations. The density of polystyrene at T = 170° C is 
p ~ 0.987 • 10 3 kg/m 3 , its viscosity is rj — 4062kg/ms, and its surface tension is 70 = 0.03 
N/m. The surface tension of Polyvinylmethylether is 70 = 0.021 N/m. A linear estimate of 
the variation of the surface tension of the mixture with the concentration gives d'y/dc = 0.018 
N/m.-S 3 - 

To calculate the coefficients of the free energy a and b,— and the coefficient of the gradient 



37 for a more detailed discussion of the 



term in the Cahn-Hilliard equation a c (cf. Ref. 
coefficients), we use the Flory-Huggins model^ 1 ^. In this model one supposes that the 
monomers of the chains of the two polymers composing the binary mixture occupy different 
sites on a square lattice. For polystyrene molecules with typically 1000 monomers, and 
a lattice spacing of ~ 10~ 10 m the random phase approximation^— gives o c ~ 10~ 16 at 
T = 170° C— This finally allows us to estimate a = —O(l) and b = 0(1). The scale of the 
concentration field C must be < 0(1) and we consider C ~ 0.5, consistent with C = \Jajb. 

The diffusion coefficients for polymers measured in experiments range from 10~ 17 to 
10~ 19 cm 2 /s.— r — We choose here the value obtained by Reiter and Steiner— of D ~ 
10~ 17 cm 2 /s because their case is close to ours. The mobility coefficient of the Cahn-Hilliard 
equation M is calculated according to the relation M = D/E with the characteristic energy 
per unit volume E = Aa 2 /b and the characteristic diffusivity D. We find M ~ 10~ 20 m 3 s/kg. 
The obtained parameters correspond to a characteristic time scale of a few hundredths of 
a second and a characteristic length scale of some nanometers. From the estimates of the 
physical parameters we obtain the dimensionless parameters Re~ 10 -17 and Ps~ 10~ 15 . 
Therefore, we set Re=Ps= but keep the ratio Re/Ps as an important parameter of or- 
der one. Normally, we will keep Re/Ps= 1 fixed but we investigate its influence below in 
Section IV CI These values correspond to the fact that polymers flow in an extremely slow 
creeping flow as expected. However, as explained before that flow has to be taken into ac- 
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count in order to explain and describe the evolution of the film profile. In consequence, the 
only free parameter that enters the bulk equations Q14p and ffl5|) is the ratio Re/Ps. Next 
we discuss the dimensionless numbers that enter the boundary conditions. For the surface 
tension number S, we obtain S~ 10 3 . However, it depends linearly on the interfacial tension 
and therefore on the type of polymers of the binary mixture. Therefore we will allow it to 
vary. Finally, a + , a - , b + , and b~ depend on the nature of the surfaces and therefore we will 
treat them as free surface parameters. 

From the original large number of parameters we, finally retain as free parameters for 
realistic binary mixtures the time scale ratio Re/Ps important for the bulk flow, and the 
numbers related to surface properties a , b^ and surface tension number S. The remaining 



dimensionless numbers discussed above or in Ref. 



371 are either very small, considered as 



zero for all purposes, or equal to one because of the scaling used. 

IV. LINEAR STABILITY OF HOMOGENEOUS FILMS 

Homogeneous base states are characterized by a concentration cq constant across the 
sample and a quiescent fluid vo = 0. Arbitrary concentrations cq are possible for neutral 
surfaces a ± = b = 0. An energetic bias at the surfaces leads to the restriction cq = 



—a + /b + = —a~/b~ (see Ref. |37J) whenever b ± ^ 0. 

Inspecting Eqs. ( 1341) to ( I4T1) we note that the evolution of the perturbation of the con- 
centration field Ci is independent of the perturbation of the velocity field u>i, i.e., the per- 
turbation fields Ci and W\ are neither coupled in the bulk equations nor in the boundary 
conditions. This decoupling of w\ and c\ results in w\ = and implies that hydrodynamics 
has no influence on the evolution of the homogeneous film in the linear stage. A further 
consequence is the lack of surface deflection defined as hi = w\j (3. 

The decoupling of velocity and concentration fields implies that all unstable modes are 
purely diffusive and correspond to the ones that can be obtained in a model based solely on 
the Cahn-Hilliard equation. We stress again that no surface deflection can develop. In the 
following we study the diffusive modes to put the present results in the context of results 
obtained in the literature. 
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A. Neutral surfaces a = b = 

Neutral surfaces correspond to the simplest possible configuration for a binary film. The 
bounding substrate and the free surface do not prefer attachment by any component. They 
do neither influence mixing or demixing phenomena. They only confine the film passively. 
Thereby, however, they restrict the possible instability modes as compared to the bulk. From 
Eq. (|34"|) we obtain the dispersion relation for neutral surfaces 

P = -q 2 {q 2 + d cc f(c)\ co ) (42) 

with q 2 = k 2 + k 2 and the vertical eigenfunctions are of the form ~ exp(ik z z). For off- 
critical mixtures (i.e., Co 7^ 0) d cc f(c)\ Co = 3cq — 1, i.e., spinodal decomposition only occurs 
for \cq\ < l/Vs. Here, we entirely focus on critical mixtures (co = 0) with 9 cc /(c)| C0 = 
— 1. As the system is infinitely extended in x-direction the horizontal wavenumber k is 
continuous. However, the vertical wavenumber can only take discrete values k z = 2irm/h 
where m = 0,1/2,1,... is a mode number. The values 0, 1/2 and 1 correspond to a 
vertically homogeneous mode, a vertical two-layer mode and a vertical three-layer sandwich 
mode, respectively. 

The m = mode is characterized by a purely horizontal variation of the concentration 
field Ci (cf. Fig. Q^a)) without any vertical structure ik z = 0). For neutral surfaces it is 
present for all film thicknesses. The corresponding dispersion relation (3 = —k 2 (k 2 — 1) is 
shown in Fig. [Ha). The maximum growth rate is (3^^ = 1/4 at A; max = 1 / -\/2. Due to mass 
conservation = for k = 0. The mode corresponds to the long-wave bulk mode of the 
Cahn-Hilliard equation.—^ 

The first mode with a vertical structure is the m — 1/2 mode. It is unstable for h > 71. 
For h < \[2tx and h > \f2ix the mode has the dispersion relations presented in Figs. Hfb) 
and[T]^c), respectively. For h < \[2tx the mode shows its maximum growth rate ff^^ 2 at 
zero horizontal wavenumber k — 0, i.e., it is purely vertical. Note that for neutral surfaces 
Pm&x^ 2 is always smaller than (equal at h — V2ir, see Fig. EJja)). For h > y/2ir the 

m = 1/2 mode has a dispersion relation shown in Fig.[]Jc) that we call 'mixed mode' because 
one has a vertical as well as a horizontal structuring (checkerboard-like) of the unstable film 
because m > and /c max 7^ 0. Note as well that /3 > at k = and /3 max = 1/4 as for 
the m = mode. As h increases the fastest growing horizontal wavenumber of the mixed 
mode becomes larger (Fig. [2](b)) and approaches the one of the horizontal mode whereas its 

14 



0.3 



CO. 



i i i i I i i i i I i i i i I i i i i I i i i i i i i i I- 




(a) 



0.2 0.4 0.6 0.8 1 
Wave number k 

Q |_[ I I I | I I I I | I I I I | I I I I | I I I I | I I I I | L 



ca 




0.1 0.2 0.3 0.4 0.5 0.6 
(b) Wave number k 

ca. ' 25 

a> 0.2 
S3 

£ 0.15 
| 0.1 
O 0.05 





(c) 



0.2 0.4 0.6 0.8 1 
Wave number k 



Figure 1: (color online) Selected dispersion relations for a homogeneous film of a critical mixture 
Co = with energetically neutral surfaces = 6 ± = and thicknesses (a) h = 3 (horizontal 
mode), (b) /i = 4 (vertical mode), and (c) h = 7 (mixed mode). The insets show sketches of the 
evolving concentration patterns for the corresponding modes. The dispersion relations discussed 
below for homogeneous films with energetically biased surfaces and as well for stratified films have 
similar shapes and will not be shown. 
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growth rate equals =1/4. The minimum of the dispersion relation at k = becomes 
deeper as h increases (thin dashed line in Fig. [2](a)). 

Higher order modes m > 1 destabilize successively at h = 2m7r. The resulting instability 
correspond to purely vertical modes with /3 max < 1/4 for h < 2^/2rmr, and to mixed modes 
with /3 max = 1/4 for h > 2\/2mix. 

Note, finally, that the dispersion relations that will be discussed below for homogeneous 
films with energetically biased surfaces and for stratified films have similar shapes as the 
ones in Figs. [T^a) to[TJ^c). However, they do not fall into the classification of modes used in 
this section. 



B. Energetically biased surfaces a^ 1 = 0, b ^ 

The relevance of energetically biased surfaces in the mixing/ demixing dynamics has been 
recognized in the design of materials based on polymer blends^. A biased surface can cause 
the formation of spinodal composition waves in a direction normal to the surfaces^ 1 ^ giving 
rise to the so-called surface directed spinodal decomposition^. 

According to the expression of the surface free energies f + {c) = 1 + a + c + b + c 2 /2 and 
f~(c) = 7 S + a~c + b~c 2 /2 and the boundary conditions ( )38|) . ( 140]) . and (14T|) . the influence 
of the surfaces enters the stability problem through the coefficients a ± that account for a 
preferential adsorption of one of the components at the respective surface, and the coefficients 

that describe mixing/ demixing properties at the surfaces that differ from the bulk. 

An energetic bias at the surfaces leads to the restriction cq = —a + /b + = —a~/b~ for 
the mean concentration.— It is very unlikely that the general case can be easily realized 
in an experiment. Therefore we restrict our analysis for the homogeneous film to the case 
6 ± 7^ with = 0. Experimental realization seems as well questionable, however, the case 
is more generic as a critical mixture yields a valid base state for any b + and b~ . The case 
was studied in depth in Ref. |67| in the context of purely diffusive decomposition of a binary 
mixture in a gap between two solid plates. This allows us to use the benchmark 
for our numerical procedure by comparing our results with the literature. 

We introduce a number of archetypal cases for b ± : the symmetric case with b + = b~ 
where both surfaces of the film have identical bias; the antisymmetric case b + = —b~, where 
the two surfaces have opposite bias; and the asymmetric case where either b + or b~ is zero. 
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Figure 2: (color online) Maximum growth rate (a) and corresponding wavenumber (b) as a function 
of the thickness for a homogeneous film of critical composition cq = and neutral surfaces = 
and = 0. Solid line: mode m = (Fig. QJa)). Dot-dashed line: mode m = 1/2 with k max = 
(Fig.QJb)). Dashed line: mode m = 1/2 with /c max 7^ (Fig. [1(c)). Thin dashed-line: growth rate 
j3 of the mixed mode m = 1/2 at the minimum of dispersion relation at /c m ; n = 0. 
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Fig. [3] presents typical results for the two most dangerous transverse instability modes 
for the symmetric (b + = b~ = 1, top row), asymmetric (b + — 1, b~ — 0, middle row), 
and antisymmetric {b + = — b~ = 1, bottom row) case. A selection of the corresponding 
eigenmodes ci(z) is given in Fig. HJa) to (c), respectively. 

In all cases one of the two most dangerous modes has a dispersion relation similar to 
Fig. QJa), i.e., /3 max — > for fc max — > 0. It is nevertheless not a purely 'horizontal mode' as 
the eigenfunction c\ shows a clear vertical structure (see Fig. H]). It dominates for symmetric 
bias (Fig. [3^a)). The other mode has normally a dispersion relation similar to Fig. [D^c), 
i-e., /3 max 7^ at A; max = 0. However, for 4.8 < h < 6 in the symmetric case /3 max occurs at 
&max = (analogous to Fig. [U(b)). This mode dominates for asymmetric and antisymmetric 
bias (Fig. [3](b) and Fig. [3|c), respectively). 

Inspecting Fig. [3] one notes that for symmetric and asymmetric bias the dependency on 
film thickness is similar to the one for the m = 1/2 mode in the case of neutral surfaces 
(Fig. [2]): All modes stabilize at a thickness h c , wavenumber and growth rate increase with 
increasing h > h c and for large h both /9 max approach 1/4. The convergence at large h is 
expected as the relative influence of the walls decreases with increasing h. Below h c the film 
is linearly stable. 

For antisymmetric surface bias (bottom row of Fig. [3]) one of the modes shows a similar 
behavior as the above modes. The other one, however, becomes progressively more unstable 
as b~ decreases (as compared to the asymmetric case). The growth rate converges for large 
thicknesses to a value much larger than the bulk value of 1/4. The corresponding eigenmode 
given in Fig. H^c) indicates that the film remains nearly homogeneous at the top where 
b + = 1 suppresses the demixing. At the bottom demixing is enforced by b~ = —1 and 
becomes much stronger than bulk demixing. 

V. RESULTS FOR STRATIFIED BASE STATES 

Non-trivial quiescent horizontally homogeneous base states of model-H are vertically non- 
uniform solutions, i.e., layered films. They are characterized by ^-dependent concentration 
profiles Co = cq(z) and vo = 0. In part I—, we discuss for different energetic biases at the 
surfaces various types of stratified film states in dependence of film thickness h. Several 
branches of stratified solutions may co-exist. They can be ordered by the mode type of 
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Figure 3: (color online) Maximum growth rate (left) and corresponding wavenumber (right) as a 
function of the thickness for an homogeneous film cq = with biased surfaces. Top row: symmetric 
bias with = 1 disfavoring demixing at both surfaces; middle row: asymmetric bias with b + = 1 
and 6~ = disfavoring demixing at the top surface; bottom row: antisymmetric bias with b + = 1 
and b~ = — 1 disfavoring (favoring) demixing at the top (bottom) surface. All solid lines correspond 
to maxima of dispersion relations similar to the one illustrated in Fig. QJa), whereas dashed lines 
normally correspond to relations similar to the one in Fig. QJc) . 



the related linear instability mode, i.e., the instability mode that destabilizes the trivial 
homogeneous solution at the bifurcation of the respective branch of stratified states. For 



details see Ref. 



37 



Although the pattern of branching is intriguing and interesting from the point of view 
of bifurcation theory, only a subset of branches has a large importance for experimental 
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Figure 4: (color online) Selection of eigenmodes c\{z) for fixed thickness ho = 5 corresponding to 
the instability modes presented in Fig. [3) (a) symmetric bias with b — 1; (b) asymmetric bias 
with b + = 1 and b~ = 0; (c) antisymmetric bias with b + = 1 and 6~ = —1. Solid lines correspond 
to maxima of dispersion relations similar to the one illustrated in Fig. QJa), whereas dashed lines 
correspond to relations similar to the one in Fig. QTc). 



systems. They can be identified calculating their energy (cf. Figs. 4, 8 and 10 of Ref. 1371 ). 
The figures show that depending on parameter values the stratified films of either type n = 
or n = 1/2 correspond to the energetic minimum. For weak surface bias, these solutions 
correspond roughly to a weak vertical stratification for n = and to strongly stratified 
two-layer configurations for n = 1/2. For stronger bias they might as well correspond to 
a 3-layer sandwich configuration. In the following, we restrict our attention to the lateral 
stability of the solutions that are energetically most favorable. 



A. Stability without hydrodynamics 

To isolate the effect of diffusive transport from the one of hydrodynamic convective motion 
on the stability of the binary mixture we first study the stability of stratified films without 
fluid motion, i.e. we set w± = in Eq. (|26|) and only solve ( 126]) . ( 13T|) . (|32|) and (|33|) together 
with the equations for the steady states. This implies as well that the free surface remains 
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Figure 5: (color online) Characteristics of the horizontal instability modes in the case of purely 
diffusive transport for neutral (a^ = = 0, solid lines in (b) and (d)) and symmetrically biased 
(a^ = 0.5, b^ = 0, dashed and dot-dashed lines in (b) and (d)) surfaces. Panels correspond to (a) 
concentration profiles cq(z) of the stratified base states at ho = 3.5 for neutral surfaces (n = 1/2 
branch, solid line) and symmetric bias (n = branch, dot-dashed line and n = 1/2 branch, dashed 
line); (c) perturbation profiles c\(z) for h = 3.5 and line styles as in panel (a); (b) maximum growth 
rate and (d) associated wavenumber as a function of the film thickness for n = and n = 1/2 base 



states (cf. Figs. 1 and 4 of Ref. 
flat. 



371 ). In all cases S= 1, Re= 0, Ps/Re= 1. 



1. Neutral surfaces 

First, we discuss the case of neutral surfaces summarized in Fig. |5] The solid line in panel 
(a) gives the concentration profile of the base state for neutral surfaces for a film thickness 
not much above the branching of the n — 1/2 branch from the trivial homogeneous state. 
The concentration decreases monotonically with vertical position corresponding to a layered 
film with component 1 enriched near the substrate and component 2 enriched near the free 
surface (remember that c = c\ — c<i). As the film is quite thin, the diffuse interface between 
the two phases nearly spans all the layer. The solid line in panel (c) gives the perturbation 
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mode ci(z). 

Figs. 13(b) and (d) show the maximal growth rate and the associated wavenumber, re- 
spectively, as a function of the film thickness. Solid lines refer to results for the energetically 
preferable base state solution. We distinguish two regimes: (i) for < h < ir the base state 
corresponds to the homogeneous film, the most dangerous mode is the horizontal mode dis- 
cussed above in Section HV Al i.e., the growth rate equals 1/4 and the horizontal wavenumber 
is constant as well (cf. Fig. [2]); and (ii) for h > n the energetically preferable base state cor- 
responds to the two-layer film shown in Fig. Efa) as solid line. It is unstable w.r.t. lateral 
perturbations, however, the growth rate and wavenumber both decay exponentially with 
increasing film thickness. The exponential decay indicates that in practical terms a film 
above 20-50 nm might appear to be stable w.r.t. lateral concentration perturbations when 
only diffusive transport is taken into account. The corresponding dispersion relations (not 
shown) are similar to Fig. Ufa), i.e., /3 max — > for & max — > 0. In the neutral case the eigen- 
modes C\(z) show an up-down symmetry (cf. solid line in Fig. E)(c)). The mean concentration 
will develop a horizontal variation as J c\(z)dz strongly deviates from zero. The variation is 
strongest along the diffuse interface. The thin solid lines in Figs. E)(b) and (d) indicate the 
behavior for the (energetically unfavorable) homogeneous n = state for h > ir. 



2. Symmetrically biased surfaces 

Results for symmetrically biased surfaces with a* = 0.5 are given as dashed and dot- 
dashed lines in Fig. [5j The energetical bias strongly influences the base states (see Figs. 3 



to 5 of Ref. l37l ). All branches known from neutral surfaces are modified already for weak 
bias. Their number reduces with increasing bias via a sequence of bifurcations. For strong 
bias normally only a 3-layer sandwich film survives. 

Even for small thicknesses h < ir the composition of the base state is not uniform any 
more, i.e., the homogeneous solutions on the n = branches are modified and become 
(weakly) stratified. The corresponding norm \\Sh\\ increases with h and a + . The strongly 
stratified n — 1/2 solution branch bifurcates before from the (modified) n = branch, 
however, the bifurcation is shifted from h = ir to h ~ 3.3 (see Fig. [5]). The bifurcated 
n = 1/2 branch is always the energetically favorable one. 

The concentration profiles given in Fig. [5](a) (dashed line for n — 1/2 branch and dot- 
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dashed line for n = branch) indicate that the preference for component 2 at the top and 
at the bottom surface causes component 1 to concentrate within the film giving rise to a 
local concentration maximum at z/h ~ 0.3 for the n = 1/2 solution and z/h = 0.5 for the 
n = solution. This corresponds to the creation of a 3-layer sandwich structure similar to 
the n = 1 solution for unbiased surfaces.— 

At small thicknesses below « it, the maximal growth rate as well as the associated 
wavenumber for the n = branch now decrease slightly with h (Fig. EJ^b) and (d)). For 
larger h growth rate and wavenumber decay faster than exponentially, until for h ~ 8.3 
the n = branch gains stability with respect to lateral perturbations. This is, however, 
practically of no importance as for h > 3.3 the n = branch is not the energetically 
favorable one. The energetically favorable one is the n — 1/2 branch that behaves similar 
as in the case of neutral surfaces. 

The eigenmodes c\ are given in Fig. |5^c). For the n = 1/2 branch the mean concentration 
c\ is smaller than in the case of neutral boundaries, i.e., the perturbation is not predominantly 
lateral. The horizontal variation is smallest at the top and strongest along the bottom 
(cf. Fig. E|c)). The n = branch is unstable w.r.t. a nearly purely lateral mode as c\(z) »s 1. 
The direct comparison of the neutral and the biased case in Fig. [5] shows, however, that the 
symmetrically biased surfaces actually decrease the lateral instability of the layered (n = 1/2) 
branch. Note finally, that the n = branch is much less unstable for symmetrical bias than 
for neutral surfaces. 



3. Antisymmetrically and asymmetrically biased surfaces 

Stability results for antisymmetric surfaces with a + = —a~ = 0.2 in the case of purely 
diffusive transport are presented in Figs. |6]^b) and (d) using dashed lines (notice that dashed 
lines are superposed to solid lines). Panel (a) and (c) present selected base state (co(z)) and 
perturbation (ci(z)) profiles. The figures give results for the n — 1/2 branch only as it is 



the energetically favorable one (cf. Figs. 6 to 8 of Ref. |37| ). 

The boundary conditions for antisymmetric surfaces favor a bi-layer structure of the 
film, i.e., they give rise to a monotonous change of cq(z) with z (Fig. |6]^a)). Growth rate and 
wavenumber of the most dangerous lateral mode both decrease with increasing film thickness 
until the film becomes laterally stable at about h = 4.5. This implies that an antisymmetric 
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Figure 6: (color online) Characteristics of the horizontal instability modes for antisymmetrically 
biased surfaces with a + = 0.2, a~ = —0.2 and = 0. Panel (a) gives base state concentration 



33) for 



profiles for solutions on the energetically favorable n = 1/2 branch (see Figs. 6 and 8 of Ref. 
selected film thickness as indicated in the legend. Panel (c) presents the corresponding perturbation 
modes c\{z) and w\{z) for h = 3.5. Note that the c\ curves with and without convection lie on 
top of each other. Panels (b) and (d) present characteristics of the horizontal instability modes 
for purely diffusive transport (dashed lines) and for coupled transport by diffusion and convection 
(solid lines). Note, that solid and dashed lines coincide for the chosen axis scales, (b) Maximum 
growth rate and (d) associated wavenumber as a function of the thickness for the branch n = 1/2. 
Parameters: S= 1, Re= and Ps/Re= 1. 

bias has a much stronger stabilizing effect than a symmetric one (compare Figs. |5]and|6]). 

The corresponding dispersion relations (not shown) are similar to Fig. [T^a). The eigen- 
modes Ci(z) for the n = 1/2 stratified film shows an up-down symmetry. It indicates that the 
concentration will develop a horizontal variation that is strongest along the diffuse interface 
and weakest at the surfaces (cf. Fig. EJ^c)). 

The stability behavior for asymmetric surface bias is shown with dashed lines in Fig. [9](b) 
and (d) for a + = 0.2 and a" = 0. The corresponding base state profiles cq{z) and eigenmodes 
ci(z) are given in Fig. [9](a) and (c), respectively. Only results for the energetically favorable 
n = 1/2 branch (cf. Figs. 9 and 10 of Ref. 1371 ) are shown. The growth rate and wavenumber 
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of the most dangerous mode decrease with increasing thickness h until the film becomes 
laterally stable at about h = 5.6. The stabilizing effect is weaker than in the antisymmetric 
case, but much stronger than in the symmetric one (compare Figs. [5] to [9]). Dispersion 
relations (not shown) are similar to Fig. [Ha). 

Before we investigate the influence of convective motion in the next section we summarize 
our findings for the purely diffusive case: For homogeneous films, we have found that ener- 
getically biased surfaces (no linear bias = 0, purely quadratic bias ^ 0) have a strongly 
stabilizing effect for small thicknesses (with the exception of an antisymmetric bias). The 
influence of the boundaries becomes weaker for larger film thicknesses. For stratified films, 
an energy bias at the surfaces (purely linear bias a ± ^ 0, no quadratic bias b ± = 0) stabilizes 
the layered film against lateral perturbations for asymmetric and antisymmetric biases. In 
both cases the n = 1/2 stratified films are stable above some critical film thickness. A sym- 
metric bias, however, does only slightly stabilize the layered films as compared to the case 
of neutral surfaces. For symmetric and neutral surfaces the maximal growth rate decreases 
exponentially with increasing film thickness, but no critical film thickness was found. 

B. Stability with hydrodynamics 

After having studied stability in the case of purely diffusive transport we now allow as 
well for transport by convection, i.e., we introduce the perturbations of the velocity fields 
back into the model. Note that the base states are identicl to the ones in the purely diffusive 
case. However, the possible convective motion of the fluid mixture may alter their lateral 
stability, and allow for an evolving deflection of the free surface, and, in consequence, lead 
to films that show lateral modulations of composition and surface profile. 

The case of a homogeneous film is discussed above in Section HVl The linear perturbations 
of the concentration and the velocity fields are decoupled implying that hydrodynamics 
has no influence on the evolution of the homogeneous film in the linear stage. No surface 
deflection can occur. 

More detailed considerations are needed for vertically stratified films as there the per- 
turbations of concentration and velocity are coupled in the linearized bulk equations and 
boundary conditions. As above we distinguish the cases of neutral and energetically biased 
surfaces. For neutral surfaces the vertical gradient of the base state concentration (<9 2 co|o,h) 
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Figure 7: (color online). Strength of the perturbation fields for neutral surfaces (a ± = 0). Panels 
(a) and (b) give the concentration field c\(x,z) with a superposed vectorial representation of the 
velocity field (vi,wi) (white arrows). Panels (c) and (d) give the corresponding stream functions 
tpi(x, z). The film thickness is (a,c) h = 3.5 and (b,d) h = 5. The color bars give the corresponding 
field 'strength'. Remaining parameters are S= 1, Re= and Ps/Re= 1. Horizontal wave numbers 
are k = 0.558 (h = 3.5) and k = 0.202 (h = 5). The lateral and vertical size of each image 
correspond to the lateral period 2ir/k and the film thickness h, respectively. 

is zero at both surfaces and increases or decreases into the bulk. This implies that the flow 
has to be driven by the internal diffuse interface between the two components [cf. momentum 
equation Eq. (JT])]. Corresponding flow and concentration perturbation fields are illustrated 
for different film thicknesses in Fig. [7J We will call this type of forcing "bulk Korteweg 
driving". 

For energetically biased surfaces the # z Co|o,/i are not zero, they might even have their 
extrema at the surfaces [cf. Fig. [5^ a)]. This implies that the driving Korteweg forces can 
be localized near the surfaces. Liquid motion, however, is suppressed at the solid substrate 
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Figure 8: (color online). Presented are the perturbation fields for symmetric surface bias (a^ = 0.5). 
Panels (a) and (b) give the concentration field c\(x, z) with a superposed vectorial representation 
of the velocity field (v\ , w\ ) (white arrows) . Panels (c) and (d) give the corresponding stream 
functions tpi(x, z). The film thickness is (a,c) h = 3.5 and (b,d) h = 5. The color bars give the 
corresponding field 'strength'. Remaining parameters are S= 1, Re= and Ps/Re= 1. Horizontal 
wave numbers are k = 0.594 (h = 3.5) and k = 0.213 (h = 5). Both solutions are on the n = l/2 b 
branch. The lateral and vertical size of each image correspond to the lateral period 2~K/k and the 
film thickness h, respectively. 

due to the no-slip condition, but can be strong at the free surface. Corresponding flow and 
concentration perturbation fields are illustrated for different film thicknesses in Fig. [SJ This 
"surface Korteweg driving" corresponds to the classical Marangoni driving at a free surface 
for one-component systems. 

Note, that strictly speaking the flow is only driven by the Korteweg term in the bulk equa- 
tions due to the exact cancellation of the Marangoni and Korteweg terms in the tangential 
stress boundary condition for the perturbations. As one is, however, still able to distinguish 
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Figure 9: (color online) Results for asymmetrically biased surfaces with a + = 0.2 and a = = 0. 
Panel (a) gives base state concentration profiles on the energetically favorable n = 1/2 branch (see 



Figs. 9 and 10 of Ref. l37l ) for selected film thickness as indicated in the legend. Panel (c) presents 
the corresponding perturbation modes c\(z) and w\. Note that the c\ curves with and without 
convection can not be distinguished by eye. Panels (b) and (d) present for for this n = 1/2 branch 
the maximal growth rate and the associated wavenumber, respectively, of the horizontal instability 
modes for purely diffusive transport (dashed lines) and for coupled transport by diffusion and 
convection (solid lines). The remaining parameters are S= 1, Re= and Ps/Re= 1. 

a driving at the free surface and a driving at the diffuse interface we call the former either 
"surface Korteweg driving" or "Marangoni driving" and the latter "bulk Korteweg driving". 
In the following we discuss the individual cases in more detail. 



1. Neutral surfaces 

The Marangoni number is expressed as Ma= a + S= a + jq/IE. Therefore an upper neutral 
surface (a + = 0) means that gradients of concentration along the free surface do not give rise 
to gradients of surface tension, i.e., do not drive a flow. Fig. [TOl shows the maximum growth 
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Figure 10: Characteristics of the horizontal instability modes in the case of coupled transport by 
diffusion and convection for neutral surfaces ( a 1 * 1 = = 0). (a) Maximum growth rate, and (b) 
associated wavenumber as a function of the film thickness for n = and n = 1/2 base states (see 



Figs. 1 and 4 of Ref. 



371 ), Remaining parameters are S= 1, Re= and Ps/Re= 1. 



rate and associated wavenumber as a function of the film thickness. A first comparison to 
Fig. [5] seems to indicates that convective transport does only cause small changes. However, 
a careful analysis shows that the instability is accelerated by up to ~ 50% (as measured at 
h = 6: purely diffusive /3 = 1.4 x 10~ 4 ; with convection (3 = 2.1 x 10 -4 ). This is difficult 
to discern in the figures as the growth rate decreases exponentially with increasing film 
thickness. 
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The effect of convective transport on the growth rate is, however, a result of a rather 
dramatic change of the velocity field inside the film. By definition there is no such field in the 
purely diffusive case, whereas now with the incorporation of convective transport 'convection 
rolls' appear driven by the bulk Korteweg term in Eq. fl27|) . This new dynamics inside the 
film allows for a surface deflection to evolve. Figs. [71^c,d) illustrate the convective motion by 
showing the stream function corresponding to the perturbation velocity field (vi,W\) in a 
vertical cut through the film for (c) h = 3.5 and (d) h = 5. In the latter case one can clearly 
appreciate how the convective cells are driven from within the bulk. At each lateral position 
there is a well visible upper cell and a less visible lower cell that rotate in opposite direction. 
For h = 3.5 the lower cell is very weak as it is strongly 'damped' by the no-slip condition at 
the substrate. The strength of the vertical flow at the upper surface is proportional to the 
evolving local surface deflection. Figs. E)^a,b) give the corresponding perturbation fields for 
the concentration C\ and as well indicate the velocity field as superposed white arrows. For 
h = 3.5 the field c\ shows nearly vertical stripes, i.e., it is almost a purely horizontal mode. 
For h = 5 the lateral modulation in the concentration field is less developed at the surfaces. 
It is strongest in the region of the diffuse interface. 



2. Symmetrically biased surfaces 

Fig. E] gives the main stability result for symmetrical bias (solid lines) with a ± = 0.5, 
and S= 1, therefore Ma= 0.5 in comparison to the case of purely diffusive transport (dashed 
lines). There are two important differences: (i) Up to h « 6 the n = branch is with 
convection only slightly more unstable. In contrast, for h > 6 the decrease in the growth 
rate slows strongly down, i.e., with further increasing thickness the branch becomes orders 
of magnitude more unstable than without convection. This is, however, only of secondary 
importance as we expect the branch to be experimentally not of great importance. It repre- 
sents the energetically favorable solution only below h « 3.3. The behavior is nevertheless 
interesting and we will later on encounter similar results for relevant branches. 

(ii) Remarkable is the splitting in two of the n — 1/2 branch. It is a consequence of 
the breaking of the up-down symmetry by the hydrodynamic boundary conditions (no- 
slip vs. free surface). Symmetries and resulting multiplicities of branches of base states are 



discussed in 



371 . In particular, for symmetric bias the n — 1/2 branch represents two solutions 
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Figure 11: (color online) Characteristics of the horizontal instability modes for purely diffusive 
transport (Dashed lines) and for coupled transport by diffusion and convection (solid lines). Surfaces 
are symmetrically biased with a ± = 0.5. Panel (a) shows the maximal growth rate and (b) the 
associated wavenumber as a function of the thickness for n = and n = 1/2 base states (see Figs. 3 



and 4 of Ref. 



371). Remaining parameters are S= 1, Re= and Ps/Re= 1. 



related by the symmetry Cq(z) — > —c^ihjl — z). The two solutions behave identically when 
increasing the energy bias. Without hydrodynamics they have as well identical stability 
properties. However, as convective transport is included the symmetry is broken and the 
two solutions acquire different stability properties (marked as branches l/2 a and l/2 b in 
Fig. [IT]) . Both branches are more unstable with convection than without, however, whereas 
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for the branch n = l/2 b the difference reaches one order of magnitude already at h ~ 5, the 
branch n = l/2 a is only slightly more unstable. 

The velocity field inside the film is quite different from the one for neutral surfaces given 
above in Fig. [7J For symmetrical bias the surface tension depends on concentration and a 
Marangoni driving at the free surface is possible. Remember, however, that the 'classical' 
solutal Marangoni effect manifests itself in the linear analysis through a surface Korteweg 
driving active in a diffuse region in the vicinity of the sharp free surface. The corresponding 
perturbations of concentration and velocity fields are illustrated in Fig. [8] for the n = l/2 b 
branch at h — 3.5 (left) and h = 5 (right). In the representation of the velocity field as 
stream function ipi(x, z) one can clearly discern 'convection rolls' driven by the free surface, 
i.e., at each lateral position there is only one convection cell that is concentrated near 
the free surface. As the driving occurs where the film is hydrodynamically most mobile the 
destabilizing effect is much stronger when caused by the Marangoni mode (surface Korteweg 
mode) than when caused by the bulk Korteweg mode. In consequence, for larger thickness 
the roll reaches less deeply into the layer. Figs. IHa,b) give the corresponding perturbation 
fields for the concentration c\ and as well indicate the velocity field as arrows. For h = 3.5 
the field c\ shows again a rather horizontal mode. It has, however, a stronger vertical element 
than in the neutral case above. Here, the lateral modulation is strongest at the substrate 
and weakest at the free surface. For h = 5 there is nearly no lateral modulation at the free 
surface. The modulation is strongest around the diffuse interface and is less developed at 
the substrate. 

As in the neutral case the convection allows for a surface deflection to evolve. Fig. [T27 a) 
shows - using h = 3.5 as example - the relative strength of the surface deflection hi = 
Wi/P as a function of the lateral wavenumber where w\ is the value at the vertical velocity 
component at the free surface. Fig. IT27 a) gives as well the corresponding dispersion relation 
j3(k) to facilitate the identification of the physically most relevant value of hi at the maximal 
growth rate. Inspecting the full range of k one finds that the maximal deflection occurs at 
small but finite k. From the maximum value hi decreases monotonically until reaching 
negative values at large enough k. Note that the mode of maximal deflection does not 
correspond to the most dangerous mode. We add two remarks regarding the interpretation of 
the hi(k) curve: (i) The fact that hi(k = 0) ^ does not contradict material conservation as 
the k = mode is marginally stable (/3 = 0) and not unstable, (ii) The deflection amplitude 
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Figure 12: (color online) Panel (a) shows the relative strength of the surface deflection hi = w\/ (3 
(solid line) and growth rate (dispersion relation, dashed line) as a function of the wavenumber for 
a fixed thickness h = 3.5 for the branch n = l/2 b . (b) The deflection h\ at the maximum of the 
dispersion relation is given as a function of the film thickness for the n = branch (solid line) and 
the branches n = l/2 a (dotted line) and n = l/2 b (dashed line) for symmetrically biased surfaces 
with = 0.5 and = 0. The remaining parameters are S= 1, Re= 0, Ps/Re= 1. 
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hi does not correspond to an absolute amplitude value, it rather represents the relative 
importance and the relative sign of the surface deflection as compared to the perturbation 
of the concentration field (cf. Section HTT|) . 

The accompanying Fig. Pl2Tb) presents the dependence of the surface deflection of the 
most dangerous linear mode on film thickness for all branches discussed above. For the 
n = branch in the range below h ps 3.3, where it is energetically favorable the deflection 
first increases from zero to ~ 0.35, then decreases again. Above h ~ 3.3, h\ continues to 
decrease till ^0.17at/i~4.2 before increasing strongly at about h ps 6 in accordance with 
the change in the f3 miLX (h) dependence discussed in connection with Fig. [TTJ For the strongly 
stratified n = 1/2 branches the deflection increases quickly with film thickness above their 
bifurcation from the n = branch indicating the importance of the surface deflection for the 
lateral instability modes. All these findings strengthen the above conclusion that at small 
thicknesses diffusion dominates even if convection is possible. Then there occurs a rather 
sharp transition over a well defined /i-range (for the present parameter set at about h — 6) 
where the convective influence becomes important. 

3. Antisymmetrically and asymmetrically biased surfaces 

Results for antisymmetrical bias with a + = —a~ = 0.2 are included above in Fig. [6] 
for the energetically preferred branch n = 1/2. Contrary to symmetrical bias, one finds 
that cases without and with convection are nearly identical. For instance, at h — 3.5 the 
growth rate with convection is only two per cent larger than without hydrodynamics and 
the ci perturbation fields can visually not be distinguished (cf. Fig. EJ^c)). In both cases 
the films stabilize w.r.t. lateral instability modes at h ~ 4.5. The following consideration 
might, however, help to understand the finding: The film stabilizes at about h = 4.5. 
Below this rather small thickness the results with and without convection are neither well 
distinguished in the other studied cases. That means the antisymmetrical bias exercises 
such a strong stabilizing influence that films are already stable 'before' (in terms of film 
thickness) convective motion becomes relevant. The velocity field W\ (cf. Fig. |6^c)) changes 
non-monotonously with z indicating a pair of convection rolls as in Fig. This indicates 
that the weak convection is driven from the diffuse interface within the film. 

Finally, we consider the case of asymmetrical bias. Fig. [9] above gives results for the 
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energetically favorable n— 1/2 branch at a + = 0.2 and a~ = in comparison to the purely 
diffusive case. The overall behavior resembles strongly the one described for symmetrical 
bias above in Section IVB2I for the influence of convection on the n = branch: up to 
h ~ 4.5 the profile is only slightly more unstable with convection than without. In contrast, 
for h > 4.5 the film is much more unstable with convection. Without convective motion 
the film stabilizes at h c ~ 5.8. Allowing for convection and therefore for the evolution 
of a surface deflection the film remains laterally unstable above h c . The growth rate and 
wavenumber of the most unstable mode still decrease roughly exponentially with increasing 
film thickness. 

The rationale behind this finding is that convective motion does not only add a second 
transport process to the dynamics but it as well allows the film to realize a different class of 
solutions, namely, films with laterally modulated concentration and thickness profiles. This 
qualitative change is as well reflected in a change of the character of the perturbation velocity 
profile (cf. Fig. EJ^c)). At h = 3.5 it shows non-monotonous behavior indicating driving from 
the diffuse interface inside the film whereas at larger thicknesses (shown for h = 5) driving 
comes from the free surface. This corresponds to the transition from bulk Korteweg driving 
with perturbation fields that resemble Fig. [7] to surface Korteweg (or Marangoni) driving 
with perturbation fields that resemble Fig. [HJ The relative strength of the surface deflection 
hi (not shown) shows similar behaviour as for the n = branch in the symmetric case (see 
solid line in Fig. [TST b)). 

We conclude this section with a remark on the character of the dispersion relation. For 
larger film thicknesses we find wavenumber ranges corresponding to complex modes. Al- 
though they seem to be of minor importance (as we find the most dangerous mode to be 
always real) they are interesting from a theoretical point of view. They seem to be closely 
related to the possibility of surface deflection, however, their physical interpretation remains 
elusive. They appear in parameter ranges where several existing real modes re-connect in a 
different way. However, we carefully checked that they are no numerical artifact: they do 
not depend on details of the used discretization and are reproducible over several orders of 
magnitude of numerical tolerances in the employed continuation procedure. Open questions 
related to this point will be the subject of further investigations. 
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Figure 13: Dependence of the maximal growth rate as a function of (a) the Reynolds number Re, 
(b) the surface number S, and (c) the ratio Re/Ps. The curves are obtained for a stratified film of 
the branch n = 1/2 with h = 3.5 and neutral surfaces = = 0. 

C. Influence of Re, Re/Ps and S 

One can expect that for different types of binary mixtures such as n-pentane/n-decane 
mixtures, polystyrene/cyclohexane mixtures, or ionic binary mixtures such as triethyl n- 
hexyl borate in diphenyl ether the parameters that we have kept fixed up to here, e.g., the 
Reynolds number Re, the ratio Re/Ps, and the dimensionless surface tension S will strongly 
differ. A full scale study of the influence of those parameters is beyond the scope of the 
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present work. However, we would like to indicate the tendencies that are to be expected by 
looking at one particular example. We show in Fig. [TBT a) to (c) the influence of changing 
Re, S and Re/Ps, respectively, on the growth rate of the most dangerous lateral instability 
mode for the case of neutral surfaces and a film thickness of h = 3.5. Wavenumbers are only 
discussed in the text. 

When changing the Reynold number (Fig. [TBTa); keeping the other parameters including 
Re/Ps fixed), the growth rate at first (for Re< 14.5) slightly decreases before increasing until 
it saturates at about Re= 1000. The corresponding wavenumber increases monotonically 
with Re and also saturates for Re ~ 1000. Note, however, that the differences between 
the minimal and maximal value of the growth rate and the wavenumber in the studied Re 
range are only about 0.3% and 0.5%, respectively. Therefore, in the present scaling and the 
studied case the variation of the Reynolds number alone has practically no influence on the 
stability of a mixture with unbiased surfaces. Remember, that the same velocity enters the 
definition of Re and the time scale that enters /3. Fixing Re/Ps furthermore implies that 
the ratio of diffusive and convective transport remains constant. 

When increasing the surface tension number (Fig. [TBT b)) the growth rate and the 
wavenumber both increase monotonically. The growth rate saturates at about S~ 5 and 
the maximum value of the growth rate and wavenumber are close to the ones found for high 
Re. Although the increase of S has a slightly larger impact than the increase of Re, the 
overall variation is still quite small: growth rate ~ 10% and wavenumber « 3%). 

The ratio Re/Ps has a much larger impact on the stability of the mixture (c.f. Fig. ITBT c)): 
The increase of the growth rate is orders of magnitude larger than when increasing Re or 
S. The growth rate increases monotonously with increasing Re/Ps, whereas the associate 
wavenumber decreases. As discussed above Re/Ps corresponds to the ratio of the typical 
velocity of the viscose flow driven by Korteweg stresses and the typical velocity of diffusive 
processes. Therefore, an increase of Re/Ps shifts the transition towards a larger influence 
of convective transport to smaller film thicknesses and results in the behavior presented in 
Fig. [13(c). 
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VI. CONCLUSIONS 



We have studied the linear stability with respect to lateral perturbations of homogeneous 
and layered films of polymer mixtures that have a free deformable surface and are bound on 
the other side by a rigid solid substrate. The paper represents the second part of a series 
of works that develops and applies a version of model-H suitable for the study of confined 
systems involving free surfaces. 

In the first part (Ref. 1371 ) we had derived a generalized model-H coupling transport 
equations for momentum, density of one component of the mixture and entropy. The model 
was then simplified for isothermal systems. Furthermore, we had modeled an energetic bias 
towards one of the components at the surfaces, had discussed how to include the evolving free 
surface and had finally determined stratified film base states for various types of energetic 
bias at the surfaces. 

The present second part has focused on the determination of the linear stability properties 



of the stratified film states studied in Ref. |37J. To this aim, we have linearized the governing 
equations with respect to small perturbations in the velocity and concentration fields, and 
the surface profile. We have emphasized that hydrodynamic flow has to be accounted for 
even in the case of extremely slow creeping flow as otherwise surface deflections can not 
evolve. The resulting linear system of equations and boundary conditions plays a similar role 
as the coupled Orr-Sommerfeld-type equations for velocity and temperature perturbations 
describing the linear surface-tension driven instability of a horizontal liquid layer with a 
deformable free-surface of a heated layer of simple liquid^. 

We have performed a numerical analysis of the stability with respect to lateral modes 
of homogeneous and stratified base states (quiescent films) for several different cases of 
energetic bias at the surfaces, corresponding to linear and quadratic solutal Marangoni 
effects. In passing we have elucidated the close relation between surface Korteweg driving 
(forcing by Korteweg stresses in the bulk momentum equation for the diffuse layer adjacent 
to the free surface) and Marangoni driving. The numerical analysis has been performed to 
a high precision relying on continuation techniques. The latter have, in particular, allowed 
to follow the most dangerous linear mode over a wide range of a number of important 
parameters. 

For neutral surfaces two types of steady film solutions exist: homogeneous and stratified 
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ones. For homogeneous films, perturbations of the composition and velocity fields are linearly 
decoupled. In consequence, the stability of the films is ruled by the linear modes of a Cahn- 
Hilliard equation in a confined slab-type geometry. We have found that the film is laterally 
unstable for all thicknesses in agreement with literature.— 1 ^ For layered films, the films 
become exponentially less unstable with increasing thickness in the cases without and with 
inclusion of convective transport. In the latter case the decrease is slightly slower, i.e., 
hydrodynamics further destabilizes the films. We have found that for neutral surfaces and 
critical mixtures the convective motion is driven by a bulk Korteweg effect, i.e., a forcing 
through the Korteweg stresses in the diffuse interface region between the two components 
of the mixture. 

Homogeneous films are still a base state for quadratically energetically biased surfaces. We 
have investigated the benchmark for comparison with literature^^^. As expected 

we have found that both - symmetrical and antisymmetrical energetic bias - stabilize the 
lateral instability below respective threshold thicknesses. This is, however, not the case for 
an asymmetrical bias; a question that needs further investigation. 

We have found that the stability behavior of layered films under different types (symmet- 
rical, antisymmetrical and asymmetrical) of energetic bias is rather rich. For purely diffusive 
transport, depending on the type of branch and type of bias an increase in film thickness 
either (i) exponentially decreases the lateral instability or (ii) it entirely stabilizes the film 
when reaching a threshold thickness. Including convective transport leads in many cases to a 
small or strong destabilization as compared to the purely diffusive case without changing the 
behavior qualitatively. Most remarkably, however, in some cases the inclusion of convective 
transport and the related widening of available configuration space for the film (it may then 
change its profile) changes the stability qualitatively. For instance, an asymmetrical energy 
bias results in a stabilization at a threshold thickness when no convection is allowed. With 
convection the film remains unstable well above this threshold, although it still becomes 
exponentially more stable with increasing film thickness. 

Finally, we have presented results regarding the dependence of the instability on several 
other parameters, namely, the Reynolds number, the Surface tension number and the ratio 
of typical velocities of convective transport driven by Korteweg stresses and by diffusive 
transport. We have identified the latter parameter as the most influential one. 

Note that we have entirely focused on the case of a critical mixture, i.e., the case of 
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zero mean concentration. A full scale analysis for off-critical systems, for instance, of the 
dependence of stability on mean concentration has been outside the scope of the present 
paper and will be pursued elsewhere. 

In general, our results have shown that the possibility of a change in the height profile of 
the fill alters the stability behavior of a film of a binary mixture not only quantitatively but 
as well qualitatively. Therefore such a possibility has always to be taken into account. 

Although, the performed stability analysis yields results regarding stability thresholds, 
and time- and length-scales of the most dangerous modes, it does not allow to predict the 
nonlinear short- and long-time evolution of an unstable film of a binary mixture. Results 
obtained on related two-layer films using long-wave sharp interface models^— show that the 
nonlinear behavior may be quite unexpected. An example are the morphological transitions 
that may occur in the process of coarsenings^. To fully account for the behavior the free 
surface model-H should be used to study the time evolution of films of binary mixtures in the 
nonlinear regime. An alternative approach could focus solely on the final stable structures, 
i.e., steady state films that are characterized by a steady non-flat surface profiles and internal 
non- homogeneous concentration profiles. 

Finally, we would like to point out that the given analysis only presents a first step 
towards a fully operative thin film model that is able to describe in a quantitatively correct 
way the coupled dewetting and decomposition observed in many experiments using films of 
polymer blends of thicknesses well below lOOnm.— . A basic ingredient of thin film physics 
- the effective molecular interactions between film and substrate - has not yet been taken 
into account. The approach followed in thin film models of simple liquids is to include 
those interactions as an additional pressure term that depends on film thickness and models 
wettability.-^. 

Nevertheless, the present theory allows to formulate a hypothesis that may explain why 
rather thick films (100-200nm) of polymer blends can rupture rather fast:— The presently 
investigated decomposition-caused convection-mediated lateral instability leads to a growing 
lateral modulation of the film surface. This allows to 'bypass' the rather slow linear phase 
of the dewetting instability as the film is taken directly into the non-linear regime of it. 
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